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We investigate differential in-plane and out-of-plane flow 
observables in heavy ion reactions at intermediate energies 
from 0.2 ~ 2 AGeV within the framework of relativistic BUU 
transport calculations. The mean field is based on microscopic 
Dirac-Brueckner-Hartree-Fock (DBHF) calculations. We ap- 
ply two different sets of DBHF predictions, those of ter Haar 
and Malfliet and more recent ones from the Tiibingen group, 
which are similar in general but differ in details. The lat- 
ter DBHF calculations exclude spurious contributions from 
the negative energy sector to the mean field which results in 
a slightly softer equation of state and a less repulsive mo- 
mentum dependence of the nucleon-nucleus potential at high 
densities and high momenta. For the application to heavy 
ion collisions in both cases non-equilibrium features of the 
phase space are taken into account on the level of the effec- 
tive interaction. The systematic comparison to experimental 
data favours the less repulsive and softer model. Relative to 
non-relativistic approaches one obtains larger values of the 
effective nucleon mass. This produces a sufficient amount of 
repulsion to describe the differential flow data reasonably well. 

I. INTRODUCTION 

Relativistic heavy-ion collisions have been extensively 
investigated to determine the nuclear equation-of-state 
(EOS) far away from saturation and at finite tem- 
perature, using semi-classical transport models of the 
Boltzmann-type The nuclear EOS which enters 

into a transport description via density and momen- 
tum dependent mean fields has mostly been based on 
phenomenological considerations by adjusting the pa- 
rameters to nuclear matter saturation properties and 
to the momentum dependence of the empirical nucleon- 
nucleus optical potential Such parametrisations, as 
e.g. Skyrme-type potentials imply different extrap- 
olations to high and low densities and high momenta, 
which should be tested in heavy ion reactions. Thus 
there have been many efforts to determine the density 
and momentum dependence of the nuclear mean field by 
studying the different aspects of the collective flow in in- 
termediate energy heav y io n collisions between 0.1 — 2 
AGeV (SIS-energies) 

On the other hand, it is well known that hadrons gen- 
erally change their properties in the medium. This basic 
feature is already incorporated in the simplest version of 
a relativistic hadronic model for nuclear matter, namely 
linear Quantum Hadron Dynamics (QHD) |]l6|, where 
the effective nucleon mass drops with density. To ob- 
tain a reasonable compressibility scalar self-interaction 
terms are introduced and finite nuclei are well described 



|l7| . At much higher densities, i.e. for the description of 
neutron stars, also non-linear terms of the vector mean 
field are required These diff'erent treatments reflect 
the inherent uncertainties in density extrapolations away 
from the saturation point. More recent and systematic 
approaches try to fix the relevant terms, e.g. by den- 
sity functional expansions of generalised QHD lagrang- 
iens which effectively incorporate the basic fea- 

tures of chiral symmetry and its breaking. However, pre- 
dictions for high densities remain questionable since effec- 
tive field theory provides a low density expansion scheme 
valid in the vicinity of the nuclear saturation point and 
below 

An alternative approach to the density and momentum 
dependence of the mean field is provided by microscopic 
many-body models. Here the nucleon-nucleon (NN) in- 
teraction is fixed by free NN-scattering and no param- 
eters are adjusted to the nuclear matter problem. In 
the relativistic Dirac-Brueckner-Hartree-Fock (DB) ap- 
proach p^-j27|] the NN-interaction is based on modern 
one-boson-exchange (OBE) potentials |Q and the in- 
medium ladder diagrams are summed self-consistently. 
This approach describes the nuclear matter saturation 
properties very reasonably, albeit not perfectly. There 
exist good arguments why DBHF should be still reli- 
able at higher densities. As pointed out in |2^ the 
nonlocal structure of the OBE potentials accounts al- 
ready effectively on the two-body level for some features 
of many-nucleon terms, e.g. intermediate A-excitations, 
and higher order effects cancel to large extent. Then the 
success of non-linear QHD lagrangiens can be understood 
in that these models parametrise phenomenologically the 
density dependence of the microscopic DB predictions 
po| . However, the constraints from finite nuclei on the 
explicit form of the different fields are limited since the 
single particle potential results from the cancellation of 
large scalar and vector potentials. Only the spin-orbit 
interaction allows to constrain the magnitude of the ef- 
fective mass On the other hand, in energetic 
heavy ion reactions scalar and vector fields are decou- 
pled by their different Lorentz transformation properties, 
in the sense that they can be tested independently, and 
additional information on the structure of the potential 
can be obtained. Thus within the DB framework there 
is the chance to attempt a unified description of different 
nuclear systems, i.e. free NN-scattering, nuclear matter, 
finite nuclei and heavy ion collision. 

However, relativistic Brueckner calculations are not 
straightforward and the approaches of various groups 
p3|-p6[ are similar but differ in detail, depending on solu- 
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tion techniques and the particular approximations made. 
In the present work we therefore want to study the DB 
predictions for the nuclear mean field in more detail in 
heavy ion collisions. In previous studies a qualita- 

tive agreement with collective flow observables has been 
found. Recently more detailed experiments have been 
performed for differential components of the collective 
flow [pT| . In particular, the rapidity and transverse mo- 
mentum dependence of collective flow has attracted much 
theoretical interest because of their strong sensitivity to 
the momentum dependence of the nuclear mean field [ p2| . 
Here we test the self energies, i.e. the nuclear EOS, 
from two different DB calculations, from ter Haar and 
Malfliet |2j] and from a more recent study performed by 
the Tiibingen group p5| , p6|l . 

In heavy ion collisions a further difficulty arises due 
to the non-equilibrium features of the phase space con- 
figurations. This has been discussed extensively in ref. 

In a fully consistent treatment, one would have 
to solve the coupled set of DB equations for the effective 
interaction in non-equilibrium nuclear matter configura- 
tions simultaneously with the kinetic equations for the 
evolution of a phase space distribution psf . However, 
such a procedure has not been realized yet, and further 
approximations are necessary in heavy ion collisions. In 
the local density approximation (LDA) the nuclear mat- 
ter mean fields are directly used in the transport calcula- 
tion. However, this approximation is not reliable enough 
at intermediate energies, because the local momentum 
space is highly anisotropic during a large part of the 
heavy ion collision ||l^,^,^. A better approximation is 
the colliding nuclear matter (CNM) approach, where the 
phase space anisotropies are parametrised by two inter- 
penetrating nuclear matter currents, i.e. the local mo- 
mentum space is given by two Fermi-spheres, or covariant 
Fermi-ellipsoids, with given Fermi momenta and a rela- 
tive velocity jsj]. In Ref. a method was developed 
to extrapolate nuclear matter DB results to CNM con- 
figurations. The CNM self energies are applied to heavy 
ion collisions in the Local (phase space) Configuration 
Approximation (LCA) ||l^ , |l5| , ^ , where the anisotropic 
phase space is locally parametrised by a CNM config- 
uration. In this paper the LCA approximation is only 
briefly discussed, details can be found in |l|,|l5|,|4) . We 
shall discuss the density and momentum dependence of 
the DB mean fields and their application to the CNM 
approximation in terms of an effective equation of state 



This paper is organised as follows: The DB formalism 
and the different DB models are reviewed in sec. II. The 
CNM approximation is discussed in sec. III. In section 
IV the transport equation and their numerical realization 
is outlined. Then we discuss the density and momentum 
dependence of the DB mean fields in ground state nu- 
clear matter, and the LCA approximation for anisotropic 
phase space configurations in heavy ion collisions. Sec- 
tion |^ contains the results of transport calculations based 
on the DB mean fields of Refs. [p5,p6| in the LCA approx- 



imation. We compare different components of collective 
fiow with respect to its energy, centrality, transverse mo- 
mentum and rapidity dependence. It is found that both 
models are able to qualitatively describe the experimen- 
tal results, but in details one finds model dependences on 
the collective flow observables. 



II. THE DB APPROACH 

Brucckner theory provides a microscopic model which 
accounts for two-body correlations in the ladder approx- 
imation in medium in the Bethe-Goldstone or the Bethe- 
Salpeter equation in the relativistic case, 



T^V + i / VGGT 



(1) 



The correlations of the Green functions, or wave func- 
tions respectively, are shifted to the effective in-medium 
interaction, i.e. the T- matrix (or CJ-matrix) [^3| . The 
in-medium propagator obeys a Dyson equation 



G = G" 



g"j:g 



(2) 



and is dressed by a self-energy E obtained in Hartree- 
Fock approximation from the T-matrix 



TG 



(3) 



The coupled set of equations (|l[]^) has to be solved self- 
consistently. In this procedure the bare interaction V, 
iterated in the Bethe-Salpeter equation, is sandwiched 
between dressed in-medium spinors. This feature is ab- 
sent in non-relativistic approaches and introduces an ad- 
ditional density dependence which is responsible for the 
signiflcantly improved saturations properties compared 
to non-relativistic tj-matrix calculations [^,^. The 
non-relativistic Brueckner approach leads to too large 
saturation densities (e.g. ps^t — 0.2 fm-^ in 13 and 
Psat — 0.24 fm^"^ in [^) and predicts a rather small 
compressibflity (K=180 MeV in HJI^). The introduc- 
tion of 3-body-forces can in principle improve on this in 
the non-relativistic case 



In contrast to the phenomenological approaches like 
QHD ^ and to effective fleld theory ||l|,|2l| this ap- 
proach is essentially parameter free. The only free pa- 
rameters are those of the NN-interaction, namely those 
of the realistic OBE potentials, which are fixed by the free 
scattering problem 1 28 1 . The success of the DB model in- 
dicates that it already accounts for the most important 
set of diagrams to describe nuclear matter. A further in- 
clusion of particle-hole correlations around the DB mean 
field also can ensure thermodynamical consistency in the 
form of the Hugenholtz- van-Hove theorem . 

In the present work we employ the results of two dif- 
ferent DB calculations, those of Ref. |E4[ and the more 
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recent ones of Ref. ||2|], denoted as DBHM and DBT, re- 
spectively, in the following. We will briefly characterise 
the difference between these two calculations. One dif- 
ference is the use of different OBE potentials; in HJ] a 
version of the Groningen potential and in |2^] the Bonn 
A potential. In both cases the same set of six non-strange 
mesons with masses below 1 GeV is used and the fit to 
the NN phase shifts is of similar quality, however, the 
actual model parameters (coupling constants and form 
factors) are different. The main difference between these 
two approaches has a more complicated origin, which is 
discussed in detail in |^,^ . The DB structure equations 
are matrix equations in spinor space. To determine 
the Dirac structure of the self energy, i.e. the scalar (Sg) 
and a vector (S'^) contributions. 



the T-matrix has to be decomposed into its Lorentz com- 
ponents, i.e. scalar, vector, tensor, etc. contributions. 
This procedure is not free from ambiguities. Due to iden- 
tical matrix elements for positive energy states pseudo- 
scalar and pseudo- vector components cannot uniquely be 
disentangled for the on-shell T-matrix. However, with a 
pseudo-scalar vertex the pion couples maximally to neg- 
ative energy states which are not included in the stan- 
dard Brueckner approach. This is inconsistent with the 
potentials used since the OBE potentials are based on 
the no-sea approximation. Hence, pseudo-scalar contri- 
butions due to the one-7r exchange lead to large and spu- 
rious contributions from negative energy states. In ||2^ it 
was shown that such spurious contributions dominate the 
momentum dependence of the nuclear self-energy, and, in 
particular, lead to an artificially strong momentum de- 
pendence inside the Fermi sea. It was further demon- 
strated in 1 25 1 that the method used in |2j] fails to cure 
this problem and in p6| a new and reliable method was 
proposed to remove those spurious contributions from the 
T-matrix. 

TABLE I. Saturation properties of nuclear matter, i.e. 
Fermi-momentum kp, saturation density psat, binding energy 
per particle E/A, effective mass m* and the compression mod- 
ulus K in the DBHF calculations of H (DBHM) and Q 
(DBT). 





kp 


psat 


E/A 


m* 


K 




[fm-i] 




[MeV] 


[MeV] 


[MeV] 


DBHM 


1.343 


0.164 


-13.6 


558 


250 


DBT 


1.39 


0.185 


-16.1 


637 


230 



The saturation properties of the two DB calculations 
are given in Table 1. It is seen there that the results of 
ter Haar and Malfliet (DBHM) |||| give a slightly bet- 
ter saturation density compared to [^ (DBT) but too 
little binding. DBT, in contrast, gives a good binding 
energy and also meets the empirical range of saturation. 
The stiffness of the EOS expressed by the compression 
moduli is similar for both approaches. A significant dif- 
ference can be observed for the magnitude of the effective 
mass but both values are consistent with the knowledge 
from finite nuclei on the strength of the spin-orbit force 
(500 MeV < m* < 700 MeV) [|o). 

We also note that the rclativistic effective mass, i.e. 
the Dirac mass m* = 7\f — given in Table 1, should be 
distinguished from the effective mass mjjp, which in non- 
relativistic approaches is used to classify the non-locality 
l4|. The latter is defined by 



(4) of the mean field |l2 



(5) 



At po 



and is approximately to^r ~ fcg = -I- n 
the two models DBHM/DBT yield the following values 
m^^{kp)/M — 0.63/0.73 which can be compared to the 
parameters used in refs. | p^ , [T3t . 

The corresponding equations-of-state are shown in 
Fig.|l|. Both models have a similar density dependence, 
however, due to the higher binding energy DBT lies gen- 
erally below DBHM. At high densities the softer char- 
acter of DBT becomes a little more pronounced. For 
comparison also the widely used Skyrme parameterisa- 
tions (soft/hard, denoted by SMD/HMD) with compres- 
sion moduh of K=20G/380 MeV |Q| are shown. It is seen 
that in particular DBT is close to the soft Skyrme pa- 
rameterisation up to about 2pQ where it starts to become 
more repulsive. 
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FIG. 1. Equation of states in the DB approaches. Solid: 
DB calculations from |2^, dashed: DB calculations from 
[ p6| . For comparison also the soft /hard (SMD/HMD) momen- 
tum dependent Skyrme parameterisations with a compression 
modulus of K=200/380 MeV Q are shown. 
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FIG. 2. Scalar (top) and vector (bottom) self energies. The 
solid lines are DB calculations from 1241 , the dashed lines from 
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FIG. 3. Momentum dependence of the scalar (top) and vec- 
tor (bottom) self ener gies at p = l/2/3po. The solid lines are 
DB calculations from [ p^ and the dashed curves from . 

Fig.^ shows the density dependence of the mean field. 
In the DBT approach the fields are generally smaller 
by about 50-100 MeV compared to DBHM. The same 
trend can be seen from Fig.^ where the momentum de- 
pendence of the scalar and vector potentials at densities 
Pq — 0.16 fm"'^, 2po and 3po is shown. In both cal- 
culations scalar and vector fields decrease with increas- 
ing momentum in a similar way. Generally, the explicit 
momentum dependence is moderate at densities p < po 
but becomes pronounced at higher densities. In ref. |2J] 



constant values of S are taken for momenta below the 
Fermi momentum whereas the results of |2^] reflect the 
full momentum dependence outside and inside the Fermi 
sea. The difference in the magnitude of the fields 
and So in the two approaches can be traced back to the 
different projection schemes discussed above. With the 
correct and complete pseudo-vector description for the 
pion contributions the fields are dominated by the a and 
u! contributions and the other mesons 7r(pseudo- vector 
couphng),p, 77, i5 give only small corrections p5| , p6t . 

Fig.^ shows the real part of the optical Schroedinger- 
equivalent nucleon potential, defined as 



in nuclear matter as a function of its laboratory energy 
-Blab — — M . Note that the definition of an op- 
tical potential is not unique in the literature, e.g. in 
refs. [12 41 1 an optical potential is defined as the differ- 
ence of the single-particle energies in the medium and in 
free space U = hP — V -|- k^. The optical potential 
defined by Eq. (B) is the Schroedinger-equivalent rela- 
tivistic potential p4| and can be covariantly defined by 



Uo^t = {kl - M^2M = {{kl + I]^)2 - m2)/2M as 



a Lorentz scalar. Even a momentum independent vec- 
tor potential So (as in the mean field approximation of 
QHD) leads to a linear energy dependence of the opti- 
cal potential (^), i.e. a momentum dependence '^y^P'P^ ■ 
The explicit momentum dependence of the DB fields falls 
asymptotically as So,s ~ {A-\-B /p) which still leads to a 
linear increase of Uo-pt at large energies. As seen in Fig. ^ 
the DB model reproduces the empirical optical potential 
Q extracted from proton-nucleus scattering for nuclear 
matter at po reasonably well up to a laboratory energy of 
about 0.6-0.8 GeV. However, it is seen that the saturation 
behaviour at large momenta cannot be reproduced by DB 
calculations. In heavy ion reactions at incident energies 
above 1 AGeV such a saturation behaviour is required 
to reproduce transverse flow observables ||ll[. Thus, DB 
mean fields start to become unrealistic around 1 AGeV. 
There exists presently no microscopic relativistic calcula- 
tion which is able to reproduce this saturation behaviour 
of the optical potential. Therefore we restrict our inves- 
tigation of transverse flow observables to an energy range 
where the DB fields can be safely applied. At higher en- 
ergies one has to rely on phenomenological approaches 
where the strength of the vector potential is artificially 
suppressed, e.g. by the introduction of additional form 
factors 0. 

One should be aware that the empirical optical po- 
tential involves densities around po and does not com- 
pletely constrain the mean fields that enter in a heavy 
ion collision, which involve large values of momentum 
and density. As a common feature relativistic DB cal- 
culation show a strong and repulsive momentum depen- 
dence also at high densities |24 2^,^l[l whereas, e.g. the 
non-relativistic C/-matrix ||3^ has a much less repulsive 
high density behaviour. 
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FIG. 4. Energy dependence of the optical potential. The 
solid lines are DB calculations from 1241 , the dashed lines from 
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at high densities. On the other hand DBT is even more 
attractive than the soft Skyrme force at low densities and 
becomes more repulsive than the soft Skyrme force only 
above 2pQ. Hence the character of the mean field ob- 
tained by the two DB model calculations show essential 
differences although the corresponding EOSs are similar. 
The test of DB fields in heavy ion collisions where high 
densities (p « 2 — 3po) and momenta greater than the 
Fermi-momentum can be reached, should allow to differ- 
entiate between the two models. 



III. THE CNM APPROXIMATION 

The DB approach discussed in the previous section de- 
scribes equilibrated nuclear matter which is characterised 
by one Fermi-sphere of a given Fermi-momentum. In 
a local density approximation (LDA) these self energies 
are directly inserted into the drift term of the RBUU 
equation. However, as discussed before, local momen- 
tum space anisotropics are a characteristic feature of 
energetic heavy ion reactions. The time scales where 
such anisotropics occur are comparable with the com- 
pression phase of the process [n5 32 33] . It was found 



that the local anisotropic momentum space can well 
be parametrised by two inter-penetrating nuclear mat- 
ter currents, i.e. by two Fermi-spheres in momentum 
space [||,|3|, or the CoUiding Nuclear Matter (CNM) 
which is schematically illustrated in Fig.^. The ap- 
plication of the relativistic DB model to CNM configura- 
tions has, however, not been realized yet, but only non- 
relativistic t^-matrix calculations have been performed for 
CNM i7|. 



P/Po 

FIG. 5. Optical potential Uopt{p,k = kp) at the Fermi 
surface as a function of density. The solid lines are DB cal- 
culations from the dashed lines from [^. For compar- 
ison also the soft/hard (SMD/HMD) momentum dependent 
Skyrme parameterisations m are shown. 





In first order the strength of the repulsion in the rela- 
tivistic case is determined by the magnitude of the vec- 
tor field. As can be seen from FigJJ the two approaches 
DBHM and DBT yield similar results at moderate den- 
sities p < po but differences become substantial at high 
densities. The generally smaller fields of DBT result in a 
less repulsive potential at high densities. This behaviour 
becomes even more clear from Fig.^j where the optical 
potentials taken at the Fermi surface, i.e. at k — kp, 
are shown as a function of the density. For compari- 
son also the potentials of the soft/hard momentum de- 
pendent Skyrme parameterisations (SMD/HMD) [|j are 
given. For DBHM the interplay between density and mo- 
mentum dependence of the optical potential is similar to 
the hard Skyrme force although the EOS is less repulsive 



FIG. 6. Schematic representation of the LDA (nuclear 
matter) and LCA (colliding nuclear matter) approximations. 

Therefore in Ref. a method was developed to ex- 
trapolate DB results to such CNM configurations. The 
CNM momentum distribution is constructed as a super- 
position of the single currents, i.e. 

/l2(k*)=/l(k*)+/2(k*)+,5/ 

= e{kp,-k'Xi) + e{kp,~k;ut^) + 5f , (?) 

where fi are the covariant momentum space distributions 
of the two NM currents and Sf — — /1/2 is a correction 
term which takes into account Pauli-blocking effects in 
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the overlap region of the two nuclear matter currents. 
The corresponding CNM self energies T,s,o{k;x) depend 
explicitely on momentum and the configuration param- 
eters X = {kp-^,kp2,Vrei} (Fermi-momenta kp-^ ^ and the 
relative velocity Vrei)- Averaging the CNM configura- 
tion, Eq. (|^), over momentum leads to mean self ener- 
gies which depend only on the parameters of the CNM 
momentum space distributions. 

The consideration of anisotropy effects in the CNM ap- 
proximation leads to non-equilibrium mean fields, which 
essentially differ from those of the equilibrium case. In 
analogy to the non-equilibrium self energies we construct 
a non-equilibrium EOS p^ , i.e. an equation of state 
which depends on the CNM parameters. In order to com- 
pare the compression energies it is useful to subtract the 
kinetic energy of the relative motion of the two nuclear 
matter currents which yields the "subtracted" binding 
energy ijj'™'^ of the system. Fig.^ shows these EOS's for 
different symmetric [kp = kp^ — kp^) CNM configura- 
tions. It is seen in both models that the effective non- 
equilibrium EOS is softer compared to the equilibrium 
EOS (solid curves for vanishing relative velocity). How- 
ever, the different density and momentum dependence 
of the two DB models leads to a different magnitude of 
this softening effect, in particular with increasing relative 
velocity. 




1.5 

Pto/Psat 

FIG. 7. Subtracted EOSs for tlie CNM approximation us- 
ing the DB results of Ref. ||^ (top) and Q (bottom). The 
Different curves give the CNM EOS's for different relative ve- 
locities {vrei = 0, 0.2, 0.4, 0.6 for solid, dotted, dashed and 
long-dashed curves, respectively). 

This softening of the effective EOS in colliding nu- 
clear matter can be understood by considering that in 
the participant zone of the reaction particles which be- 
long to projectile and target are separated in momen- 
tum space. This geometrical effect works in a simi- 
lar way as an additional degree of freedom and leads 
to a softening of the effective EOS experienced by the 
nucleons in such configurations. This type of phase 
space effect is not included in standard transport cal- 



culations for heavy ion collisions, even when momentum 
dependent interactions are used Phenomenolog- 
ical mean fields C/(x, k) — J7ioc(e(x)) + C/nonioc(x, k) are 
usually composed by a local, density dependent poten- 
tial [/ioc(f?(x)) and a non-local momentum dependent 
part [/nonloc(x,k) = / ^^^'/(x, k')F(x, k - k') with V 
an effective momentum dependent two-body-interaction. 
Whereas C/nonioc(x, k) accounts properly for the actual 
momentum space configurations /(x, k'), the local part 
does not depend on the momentum space. Consequently, 
U\oc{q) reflects a density dependence which is correct 
in equilibrated nuclear matter but does not apply to 
anisotropic momentum space configurations. 



IV. THE TRANSPORT MODEL 

In the present work heavy ion collisions are treated by 
the relativistic (R)BUU equation ||,|l|,|l5| 



{k^d^^ + [klF^'' + m*(9>*)] d^; } f{a 



Icoll 



(8) 



This equation describes the evolution of a classical one- 
body phase space distribution f{x,k*) under the influ- 
ence of a self-consistent mean field, or the scalar and 
vector self energies and S'^. The self-energies de- 
termine effective momenta and masses of the dressed 
quasi-particles in the nuclear medium fc*'^ = k^ — S'^, 
TO* = M — . The field strength tensor of the vector 
field F^'^ = d^Yj'^ — d'^T,^ gives rise to a Lorentz force 
as in electrodynamics. This introduces in a most natural 
way a first order momentum dependence which in non- 
relativistic treatments has to be parameterised explicitely 
p],ppC|]. The collision term describes 2-body collisions 
and is treated by cascade-like Monte-Carlo simulations, 
as in relativistic versions of the QMD-model . We in- 
clude the relevant nucleonic excitations at SIS energies, 
i.e. the A(1232) and N*{1A40) resonances and their de- 
cay to one- and two-pion states. The cross-sections for 
elastic and inelastic scattering as well as differential cross 
sections are taken from Ref. Q| which are used also 
in QMD calculations at SIS energies |44p5|]. The drift 
term of the RBUU equation (||) is numerically treated in 
the relativistic Landau- Vlasov method (RLV) . This 
is a test particle method, where the test particles are 
represented by manifestly covariant Gaussians in phase 
space. In order to reduce numerical fluctuations a num- 
ber of 50 — 100 test particles per nucleon was found to 
be sufRcient here. The nuclei are initialised to fit density 
profiles obtained from self-consistent Thomas-Fermi cal- 
culations [|6|. For both models (DBT and DBHM) we 
adjust the testparticle distributions to the same reference 
distribution. The corresponding initialisations are stable 
over the durations of the considered reactions. Energy- 
momentum conservation is fulfilled with an accuracy of 
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3 — 5 % of the initial kinetic centre-of-mass energy of the 
colhding nuclei. 

The anisotropic phase space effects discussed above are 
incorporated in heavy ion collisions by applying the CNM 
or non-equilibrium DB mean fields in the framework of a 
Local Configuration Approximation (LCA) 33 1. In 

this approach the phase space is parametrised locally by a 
CNM configuration where the invariant configuration pa- 
rameters X directly determined from the phase space 
distribution f{x, k*). Transport calculations have shown 
that the collective flow is reduced if the non-equilibrium 
effects are taken into account which is consistent with 
a softening of the effective EOS in heavy ion collisions as 
discussed above. 



V. COLLECTIVE FLOW EFFECTS 

In this section different types of collective flow observ- 
ables are investigated in transport calculations using the 
DBHM/DBT mean fields in the local phase space config- 
uration approximation (LCA). To compare with experi- 
ments the same methods are used for centrality selection 
and reaction plane determination, and the theoretical re- 
sults are subjected to filter routines to simulate the ex- 
perimental detector efficiencies, if necessary. Since we 
mainly compare to results from the FOPI Collaboration 
|,|l]j4|0 these correspond to the FOPI (Phase-I or 
Phase-II) set up. These filters are sensitive to fragment 
distributions. Thus we also generate fragments in the fi- 
nal state of the reaction (at « 100 — 200 fm/c depending 
on the incident energy) using a phenomenological phase 
space coalescence model. The coalescence parameters in 
coordinate and momentum space are separately adjusted 
for both models DBT and DBHM to fit the experimen- 
tal mass distributions of light fragments {Z < 3) Jl5| , p0t . 
After the cluster formation we apply the same procedure 
for the reaction plane reconstruction as the FOPI group. 
The corresponding corrections are very close to those ob- 
tained with the IQMD model ||5l[ . 
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FIG. 8. ERAT cross sections for protons in Au-|-Au reac- 
tions at 0.4 AGeV beam energy. Transport calculations with 
the DB mean fields from Refs. H (DBHM) |||| (DBT) are 
compared with the FOPI data p9| . 

Criteria for the determination of the centrality class of 



an event are the multiplicity of charged particles (PMUL) 
and/or the ratio of transversal to longitudinal energy 
{ERAT) As an example, in Figs. | and | these 

observables are shown for protons at an energy of 0.4 
AGeV and compared to FOPI data ^SE^. Fig. dis- 



plays the differential cross section da/d{ERAT). Large 
ERAT values correspond to central reactions whereas 
small values indicate semi-peripheral and peripheral re- 
actions [|l5|,^,|2j. A qualitative agreement with experi- 
ment is achieved for both mean fields, with DBHM some- 
what overestimating the data. At small ERAT values 
below about 0.2 the quantity da/dERAT is strongly af- 
fected by trigger effects and drops very rapidly. As dis- 
cussed in model calculations tend to strongly over- 
estimate the data below 0.3, which is due to detector 
cuts which remove the projectile remnant. Since such 
remnants, i.e. very heavy clusters, are not formed in 
transport calculations the cross section is overpredicted 
in very peripheral collisions although the total reaction 
cross section coincides with the experimental one. 




50 

PMUL 

FIG. 9. Multiplicity of charged particles. The solid lines 
are DB calculations from the dashed lines from p^ . 

In Fig. ^ the multiplicity distributions da/d{PMU L) 
are shown. Also here the agreement with experiment is 
reasonable for DBHM and quite good for DBT. Particu- 
larly, we can reproduce the plateau (in logarithmic scale) 
of these distributions which is used in the centrality se- 
lection. Altogether, the good agreement with experiment 
makes the centrality selection reliable using either the 
PMUL or ERAT observable. In the the following we 
will mainly use the PMUL criterion. The correlation 
between multiplicity and centrality intervals is defined 
as in experiments in the following way: the lower limit 
of the highest multiplicity bin, called PM5, is fixed at 
half of the plateau value, and the remaining multiplicity 
range is divided into four equally spaced intervals, de- 
noted by PM4 to PMl. PM5 then corresponds to most 
central reactions, and PMA and PM3 to semi-central 
and peripheral ones, respectively [E^. 



A. Nuclear stopping 

We start the flow analysis with the longitudinal distri- 
butions which are characterised in terms of the rapidity 
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Ycra = iln(l 

rapidity Y^°^ -- 



— /3c. m.)- Here the normalised 
./FjP™j is considered. 
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FIG. 10. Rapidity distributions for Au+Au collisions at 
0.4 AGeV beam energy. The different centrality intervals run 
from peripheral (upper left) to central (lower right) collisions. 
The solid lines are DB calculations from the dashed lines 
from 0. 

The rapidity distributions are shown in Fig.|l^ for dif- 
ferent centraUty classes (using the ERAT selection ) . 
This observable is strongly affected by detector cuts 
p9| , ^ , which is reflected in the asymmetry of these dis- 
tributions relative to the cm. -rapidity, due to the angular 
limitations of the FOPTdetector (Phase-1) H]. 

It is seen that both models are able to generally re- 
produce the experimental results. Stopping is influenced 
most strongly by the choice of the NN cross sections, 
which are the same in both models. An indirect influ- 
ence also from the EOS can be expected since a softer 
and less repulsive EOS leads to more compression and 
thus to more collisions. However, the differences are not 
very pronounced which reflects the fact that the EOSs 
are still similar in the explored density regime and do 
not differ too much even at the maximal densities reach 
in central reactions (p « 2 — 3 po)- In peripheral col- 
lisions, where the deflection of spectator matter by the 
repulsive momentum dependent component of the mean 
field plays a more dominant role, the rapidity distribu- 
tions show bigger differences. Similar trends have been 
observed in in Ref. |^ with soft/hard Skyrme forces. 



B. In-plane flow 

Next we consider the emission of matter projected onto 
the reaction plane described by the mean in-plane or side- 
ward flow 1^,^. Fig. |ll] compares the mean in-plane 
proton flow < Px/A > in semi-central (PM4) Au-|-Au 
collisions at incident energies of 0.25, 0.4 and 0.6 AGeV 
to the FOPI data from ||5|]. Both models reproduce the 
energy dependent increase of the proton flow and are gen- 
erally in good agreement with the data. We observe that 
DBHM leads to a stronger in-plane flow in the specta- 
tor rapidity region whereas the slope near mid-rapidity 
~ is less affected by the differences in density and 
momentum dependence of the mean fields. At 0.6 AGeV 



DBHM starts to overpredict the slope of the in-plane flow 
whereas DBT is still in reasonable agreement with exper- 
iment. This reflects again the more repulsive character 
of the DBHM forces which becomes more pronounced 
with increasing energy, in particular in the spectator re- 
gion y*^"^ ^ 1. On the other hand, DBT describes the 
in-plane flow well around midrapidity but slightly under 
predicts it in the spectator region. 
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FIG. 11. Mean in-plane transverse flow of protons versus 
normalised rapidity for semi-central (PM4) Au-(-Au collisions 
at 0.25, 0.4 and 0.6 AGeV beam energy. The solid lines are 
DB calculations from the dashed lines from |^^. The 
experimental results are taken from js^ . 

The situation is more clearly seen when the mean di- 
rected flow P^^^ , i.e. the in-plane flow integrated over the 
forward hemisphere (Y^°^ > 0), is considered. In Fig.p^ 
the excitation function of P^"" is compared to the FOPI 
data ll^, again for PM4. P^"" is a measure for the over- 
all repulsion experienced by the reaction. The compar- 
ison with data shows that both microscopic models can 
explain the experimental results relatively well. How- 
ever, above 0.4 AGeV the experimental excitation func- 
tion of P*"' shows a saturation behaviour which is not 
completely reproduced. Here, as already seen in FigjTl|, 
DBT is in good agreement with experiment at higher en- 
ergies but underpredicts the flow at lower energies. For 
DBHM the situation is just opposite. Due to the FOPI 
acceptance P^*'' is dominated by the bounce-off of spec- 
tator matter at V*^°^ ~ 1 (see the corresponding rapidity 
distributions in Figjl^) which explains the low value of 
DBT at 0.4 AGeV. The more repulsive character of the 
DBHM mean flelds, in particular at high densities, pro- 
duces a bounce-off of the spectator remnants in the re- 
action plane, resulting in too high transverse momenta 
near spectator rapidities at energies above 0.6 AGeV. 
Since the maximal densities reached in the model calcu- 
lations are changing very little between 0.4 and 0.8 GeV 
the slope of the P**" excitation function is determined 
mostly by the momentum dependence of the mean field 
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and less by the density dependence of the EOS. This may 
be considered as a way to test the average momentum 
dependence of the models. However, as discussd below 
an integrated observable reflects the average dynamics 
which is dominated by densities around po- To extract 
more decisive information on the momentum dependence 
at supranormal densities differential observables should 
be considered. 
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tides experience their most violent changes in momen- 
tum. This can be done by the following quantity 
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FIG. 12. Excitation function of the mean directed in-plane 
flow for semi-central (PM4) Au-|-Au collisions (data from 
[p3|). The solid lines are DB calculations from [Q, the dashed 
lines from p^ ]. 

As a differential observable the dependence of the in- 
plane flow on transverse momentum pt has recently at- 
tracted great interest. While the global behavior of the 
transverse flow, as expressed by P^"" or < p^/A >, gives 
an average over the entire evolution of the collision, the 
Pt dependence allows to obtain information on different 
stages of the reaction. High pt nucleons, but also pi- 
ous ]4^ ] , originate from the early and high density phase 
of the reaction. This is a general feature of heavy ion 
collisions and seems to hold at bombarding energy from 
SIS @ to AGS and SPS energies [1^. In Ref. (l| is 
was pointed out that the pt dependence of the transverse 
flow in peripheral reactions is particularly sensitive to 
the momentum dependence of the nuclear mean field at 
supranormal densities. Consistent with this observation, 
in JTst it was demonstrated that the emission time of 
high Pt particles at midrapidity coincides well with the 
high density phase of the reactions. To be more quan- 
titative, in Fig. 13 we consider the correlation between 
the averaged density < ps > at which the particles are 
emitted and the transverse momentumpt . Exemplarily a 
semi-peripheral (b=6 fm) Au+Au reaction at 0.4 AGeV 
is considered. Since particles interact with the surround- 
ing matter over the entire duration of the reaction more 
or less strongly - with the mean field and by binary col- 
lisions - it is a priori not clear how to define an emission 
time and, correspondingly, a density from which the par- 
ticles carry information. It is, however, natural to select 
that time and that corresponding density where the par- 



(9) 



where the index i runs over all baryons. Eq.(^ samples 
over the densities weighted by the relative changes of the 
momenta. This is a generalisation of the "freeze-out" 
density because it takes the whole reaction history into 
account. Since experimentally only the final state mo- 
menta are detected, this quantity is useful to establish 
a correlation between final observables and the stages 
of the reaction on which the particles carry informa- 
tion. The resulting emission densities are generally mod- 
erate because one averages over the entire phase space, 
i.e. over participant and spectator regions, and over the 
complete reaction history including the dilute expansion 
phase. Applying different pt-cuts (p'f^ = Pt/Pc^^.) one 
sees that the low pt-particles and hence the bulk of the 
nucleons experiences a dynamical evolution which takes 
mainly place at densities below saturation density when 
integrated over the complete space-time of the reaction. 
The high p(-particles {p[^^ > 0,5), in contrast, are gov- 
erned by larger densities. As a consequence, they probe 
the momentum dependence of the optical potential at 
supranormal densities. In Fig. 13 as in the calcualtions 
below (Fig.|l4|) we considered projectile rapidities, but 
the observed behaviour is even more pronounced in the 
midrapidity region. 
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FIG. 13. The density distribution which is correlated with 
the dynamical history of the particles is shown for different 
Pj"' cuts for a semi-central (b=6 fm) Au-|-Au collisions at 0.4 
AGeV. 

The flow and its p^^^ -dependence has been discussed in 
terms of a Fourier analysis of the experimental azimuthal 
distribution fsHl 
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— (pf', Y^"^) = vo{l + 2vicos{cb) + 2v2Cos{2cb)) . 
acp 

(10) 

In Eq. ( p^ vo is a normalisation constant, vi{pf \ F*^"^) 

describes the in-plane collective flow and V2ip[^\ Y^^'') 
the emission perpendicular to the reaction plane, also 
called elliptic flow. The quantities vi^2 can be determined 
directly as vi =< p^/pt > and V2 =< (pi -pl)/pl >, 

where pt — '\JPx + Py the transverse momentum per 
nucleon. 

Fig. |l^ shows the -dependence of the in-plane flow 
{vi) at normalised rapidities 0.5 < 1^*^°^ < 0.7 for semi- 
central (PM4) Au+Au reactions at 0.4 AGeV incident 
energy for protons {Z = 1) and for light fragments {Z = 
2 + 3). 



comparison to data favours a weaker repulsion at higher 
densities as predicted by DBT. In contrast to the global 

Pf observable the quantity vi{p[°\ Y^°^) IS more sen- 
sitive to high density matter in the participant region. 

The situation is similar for light fragments {Z = 2 + 3) 
which, however, show a stronger collectivity compared 
to free nucleons, as also seen by the higher value of V2 
in Fig.^. This can be understood by assuming that the 
heavier fragments are mainly formed in the spectator re- 
gions, whereas the free nucleons originate from the en- 
tire phase space. Thus, one obtains on the average a 
higher mean transverse momentum < Px/A > for frag- 
ments than for nucleons. This scenario of the fragment 
formation is consistent with the findings in Ref. pQ] that 
the spectator enters into an instability region with con- 
ditions which are near the experimental ones for a liquid- 
gas phase transition fs^ . 
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FIG. 14. In- pla ne flow in terms of the first Fourier coeffi- 
cient (see Eq. (p_0|)) for semi-central (PM4) Au-|-Au collisions 
at 0.4 AGeV for protons and light fragments (data are from 
[^). Statistical uncertainties of the calculations are indi- 
cated by bands. 

A significant dependence of the differential in-plane 
flow on the DB model is observed. At low transverse 
momenta {p^^^ < 0.2) the two models do not show large 
differences. Since the mean transverse flow is dominated 
by low Pt particles the < Px/A > observable therefore 
does not differentiate very much between the two models 
in the present situation. Consequently, in Fig.|ll| the two 
curves for < Px/A > (0.4 AGeV) are very similar to each 
other in the considered rapidity range. However, the dif- 
ference between the models becomes more pronounced 
when higher transverse momenta {p^^^^ > 0.2) are stud- 
ied. In Fig.|lJ DBT reproduces the data very well for 
both, protons and light fragments, whereas DBHM over- 
predicts the flow significantly which reflects the its more 
repulsive character at baryon densities above po- The 



C. Out-of- plane flows 

A preferential out-of-plane emission of particles, the 
so-called squeeze-out, is characterised by negative values 
of the second Fourier coefficient V2 in Eq. (|lC)|) w hereas 
a positive value of V2 indicates in-plane flow ]55|]. A re- 
lated variable to characterise the azimuthal anisotropy of 
par ticle emission is the squeeze-out ratio Rn defined by 



Rn ^ {N{<i> = 90°) + iV($ = 270°))/(iV($ 0°) 

+ 7V($ 180°)) = (1 - 2t>2)/(l + 2t-2) . (11) 

In terms oi Rn an isotropic emission corresponds to 
Rn = 1, i?Ar > 1 indicates squeeze-out and Rn < 1 a 
preferential in-plane emission. At lower energies squeeze- 
out is mainly due to shadowing of the participant parti- 
cles by the spectators and it is therefore most pronounced 
for mid-rapidity particles. A sensitive way to probe the 
momentum dependence of the mean field is the transverse 
momentum dependence of the elliptic fiow V2. In contrast 
to central reactions where the squeeze-out is influenced 
by both, the static and the momentum dependent part 
of the mean fleld, peripheral reactions allow to decou- 
ple the momentum dependence to some extent from the 
static part. In Ref. |13] this has been demonstrated for V2 
at high Pt in peripheral reactions and was also confirmed 
by the investigations of ||l^. Since high pt particles al- 
most exclusively stem from the high density phase they 
probe the momentum dependence of the optical potential 
at supranormal densities, as was seen also in Fig. 13. 

First in Fig. |l^ we consider semi-central (PM4) re- 
actions. It shows the azimuthal distributions at mid- 
rapidity for all charged particles (nucleons plus frag- 
ments) in Au-f An at 0.6 AGeV. The calculations are 
compared to experimental results from FOPI |^l[] for 
Z = 1 and Z = 2 fragments. The data show a stronger 
out-of-plane emission for Z = 2 fragments relative to 
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those for nucleons which is due to to the higher collec- 
tivity of fragments as discussed above. Due to limited 
statistics the theoretical values in Fig. 15 are given for 
all charged particles. We observe a qualitatively good 
agreement between theory and experiments, although 
both calculations slightly overestimate the experimental 
squeeze-out signal. 
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FIG. 15. Azimuthal distributions at 

mid-rapidity {\AY^°^\ < 0.15) for semi-central Au+Au col- 
lisions at 0.6 AGeV beam energy (data from [E^). 
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FIG. 16. Transverse momentum dependence of the elliptic 
flow at mid-rapidity (lAF^^'l < 0.1) V2 for peripheral (PM3) 
Au+Au collisions at 0.4 AGeV. Data are taken from wA. 

Fig.|l^ shows the pt dependence of the elliptic flow V2 
in peripheral Au-f Au collisions (PM3 multiplicity inter- 
val, 5mcan ~ 7fm) at 0.4 AGeV. Consistent with the pic- 
ture that squeeze-out is mainly due to shadowing, which 
is most effective in the early high density phase of the 
reaction, the elliptic flow becomes increasingly negative 
with increasing pt- A physical explanation for this be- 
haviour is that a strong repulsive momentum dependent 
component of the nuclear interaction deflects these par- 
ticles from the central zone more violently. This leads 



to enhanced shadowing by the spectator remnants inside 
the reaction plane and a stronger emission out-of-plane. 
This effect is clearly seen for the two types of interactions 
used here: Due to the stronger repulsive momentum de- 
pendence of the DBHM forces at supranormal densities 
the squeeze-out signal increases much stronger with pt 
than for DBT, which on the other hand, is closer to the 
data of 

This behaviour is also consistent with the findings of 
Ref. where soft and hard Skyrme forces within the 
framework of the QMD model were subjected to a com- 
parison of the same observable. The softer equations- 
of-state (soft Skyrme in and DBT in the present 
case) yield a slower increase of V2 and are in reasonable 
agreement with experiment. However, in Ref. |4^ ] both 
version of Skyrme forces have an identical momentum 
dependence whereas in the present case DBHM is signif- 
icantly more repulsive at large densities. In the analysis 
of Ref. on the other hand, parameterisations were 
used which show the same density dependence, i.e. the 
same EOS, but differed in their momentum dependence. 
A stronger repulsive character of the model, expressed by 
both, a stiffer EOS and/or a stronger momentum depen- 
dence generally results in a stronger squeeze-out signal 
for high Pt particles. The authors of Ref. |jl^ favour a 
parameterisation (HM in Ref. jl^) which yields results 
for this observable close to DBHM in our case. The com- 
parison to the data of ji^ in Fig.|l^ supports, however, a 
weaker momentum dependence at supranormal densities. 
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FIG. 17. Transverse momentum dependence of 

the first {vi, left) and second (u2, right) azimuthal 
flow Fourier-coeflicients for the system Ru+Ru at 0.4 
AGeV incident energy extracted in rapidity intervals of 
-0.7 < Y'"' < -0.5 (left) and -0.3 < F*") < -0.1 (right) for 
semi-central (PM4) and peripheral (PM3) reactions, respec- 
tively (data from ]3l|). 

In Figjl^ the same analysis is performed for the Ru + 
Ru system. Here the differential components of the 
in-plane and out-of plane flow in terms of the Fourier- 
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coefficients vi and V2 as a function of the transverse mo- 
mentum and rapidity are shown. The rapidity intervals 
are -0.7 < y(°) < -0.5 and -0.3 < Y^"^ < -0.1 for 
the vi and V2 analysis, respectively. Both, DBHM and 
DBT are in fair agreement with the FOPI data |^l| on 
the V2 observable. However, the pt dependence of the 
in-plane flow vi strongly depends on the model. As al- 
ready observed for the Au + Au system (Fig.Q) DBHM 
strongly overpredicts the flow vi at intermediate pt ■ Thus 
we conclude that the differences observed in the vi ob- 
servable are due to the different momentum dependence 
of the fields. The V2 observable shows here much less 
model dependence than e.g. in the Au + Au system. In 
the smaller system the compressional effects are smaller 
which reduces the sensitivity to the different momentum 
dependence. 

In Fig. [ij the centrality dependence of the squeeze- 
out ratio Rm at midrapidity (jAF^"^! < 0.15) for parti- 
cles with high transverse momenta (0.4 < P(°^ < 0.55) 
is shown. The decrease of the squeeze-out ratio i?Ar at 
higher impact parameters can be explained by the strong 
vector repulsion of the nuclear mean field. Again both 
models can reproduce the centrality dependence of the 
experimental data qualitatively, but not in detail. The 
DBHM calculations overpredict the data at low impact 
parameters, whereas DBT underpredicts the data for pe- 
ripheral collisions. Similar results have been found in 
recent studies of the FOPI collaboration in comparison 
with the IQMD model with soft/hard Skyrme forces 
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FIG. 18. Centrality dependence of the squeeze-out ratio 
Rn at mid-rapidity (|AF'°'| < 0.15) for Au+Au collisions at 
0.4 AGeV (data from |^]). Statistical errors of the calcula- 
tions are indicated by bands. 

Finally, Fig. ^ shows the excitation function of the 
total elliptic flow V2 at mid-rapidity. With increasing 
beam energy the emission is preferentially perpendicular 
to the reaction plane {v2 < 0) which approaches a max- 
imum around 0.4 AGeV and then decreases again. This 



behaviour can be understood by shadowing and compres- 
sion effects as discussed in detail in Ref. . Again both 
models reproduce the general trend of a sample of exper- 
imental data from FOPI EOS HL Plastic BaU [|, 
LAND [| and E895 ||| shown in Fig.|l9|. However, with 
DBHM the maximum of negative V2 values is slightly 
shifted to higher energies and the absolute values are 
larger. The experimental data show strong variations in 
the magnitude of the V2 coefficient which cover the range 
of the theoretical calculations. However, DBT better re- 
produces the maximum at the correct energy and lies 
closer in absolute magnitude to the most recent data from 
FOPI Assuming that the latter are the most reli- 

able measurements, e.g. with respect to reaction plane 
corrections etc., the comparison to experiment favours 
the DBT mean fields. 
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FIG. 19. Energy dependence of the elliptic flow V2 at 
mid-rapidity. The data are taken from the FOPI- (triangles) 
[0, EOS (diamonds) H, Plastic Ball (triangles down) Q, 
LAND (triangle right) g and E895 (circle) M. 



VI. SUMMARY 

We investigated the collective nucleon flow in heavy 
ion collisions at intermediate energies (0.15 1 AGeV) 
within a relativistic BUU transport model with mean 
fields based on relativistic Dirac-Brueckner-Hartree-Fock 
(DB) calculations for nuclear matter. The anisotropy of 
the local momentum space in the participant region of 
heavy ion reactions was taken into account in the collid- 
ing nuclear matter approximation. 

We compared two different DB fields, those of ter Haar 
and Malfliet (DBHM) and the more recent calculations 
performed by the Tubingen group (DBT). From a the- 
oretical point of view the latter ones have a stronger 
physical foundation since spurious contributions from 
a strong coupling of a pseudo-scalar ttNN vertex to 
negative energy states were removed by an improved 
projection scheme for the in-medium T-matrix. The 
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DBT model EOS is slightly softer (K=250/230 MeV 
for DBHM/DBT) and the effective mass is significantly 
larger (m* = 558/637 MeV for DBHM/DBT) than for 
the DBHM model. This results in smaller fields and a 
less repulsive optical potential for DBT. The smaller re- 
pulsion of the DBT model is also expressed in a terms of 
a larger non-local effective mass (tojv_r — 0.63/0.73 M for 
DBHM/DBT) which is commonly used in non-relativistic 
approaches in order to classify the strength of the mo- 
mentum dependence of the potential. 

Both models yield a reasonable description of in-plane 
and out-of-plane flow observables. A more detailed com- 
parison to data, in particular the transverse momentum 
dependence of vi and V2 favours the softer EOS and 
the less repulsive character of the DBT optical poten- 
tial at supranormal densities. This observation is con- 
sistent with the information obtained from subthreshold 
K'^ production where also the scenario of a soft EOS 
is supported At 0.6 and 0.8 AGeV DBT also 

yields a very accurate description of the transverse in- 
plane flow whereas at lower energies the in-plane flow 
requires some more repulsion as provided by the DBHM 
model. Interestingly, a similar observation was made in 
the non-relativistic approach of Ref. . 

In summary, the microscopic DB approach, where no 
parameters are adjusted to the nuclear matter saturation 
properties nor to the empirical optical nucleon-nucleus 
potential, predicts a density and momentum dependence 
of the mean field which to a large extent is consistent 
with the observations from heavy ion collisions. 



ACKNOWLEDGMENTS 

We are grateful to A. Andronic from FOPI for many 
discussions concerning the interpretation of data. 



[1] E.A. Uehling, G.E. Uhlenbeck, Phys. Rev. 43 (1933) 552; 
P. Danielewicz, Ann. Phys. 152 (1984) 239; G.F. Bertsch, 
S. Das Gupta, Phys. Rep. 160 (1988) 189. 

[2] B. Blattel, V. Koch, U. Mosel, Rep. Prog. Phys. 56 
(1993) 1. 

[3] S. Kama, B.C. Clark, E.D. Cooper, H.S. Sherif, R.L. Mer- 
cer, Phys, Rev. C41 (1990) 2737; E.D. Cooper, S. Hama, 
B.C. Clark, R.L. Mercer, Phys. Rev. C47 (1993) 297. 

[4] J. Aichelin, H. Stocker, Phys. Lett. B176 (1986) 14. 

[5] W. Reisdorf and E.G. Ritter, Annu. Rev. Nucl. Part. Sci. 
47 (1997) 663. 

[6] N. Herrmann, J. P. Wessels, T. Wienold, Annu. Rev. 
Nucl. Part. Sci. 49 (1999) 581, and refs. therein. 

[7] T. Maruyama, W. Cassing, U. Mosel, S. Teis and K. We- 
ber, Nucl. Phys. A573 (1994) 653. 



[9: 

[lo; 

[11 

[12; 

[13 
[14 

[is: 

[16 
[17 

[is: 

[19 
[20 

[21 

[22: 
[23: 

[24' 
[25^ 

[26: 

[27 

[28: 



[29 
[30 



[31 

[32: 
[33: 



[34' 
[35: 



[36] 



Y.M. Zheng, CM. Ko, Bao-An Li and Bin Zhang, Phys. 
Rev. Lett. ASS (1999) 2534. 

J. Nemeth, G. Papp, H. Feldmeier, Nucl. Phys. A647 
(1999) 107. 

P.K. Sahu, A. Hombach, W. Cassing, M. Effenberger, 
U. Mosel, Nucl. Phys. A640 (1998) 493; A. Hombach, 
W. Cassing, S. Teis, U. Mosel, Eur. Phys. J. A5 (1999) 
157. 

P.K. Sahu, W. Cassing, U. Mosel, A. Ohnishi, Nucl. 
Phys. A672 (2000) 376. 

P. Danielewicz, Nucl. Phys. A673 (2000) 375. 

A. B. Larionov, W. Cassing, C. Greiner, U. Mosel, Phys. 
Rev. C62 (2000) 064611. 

C. Fuchs, T. Gaitanos, H.H. Wolter, Phys. Lett. B381 

(1996) 23. 

T. Gaitanos, C. Fuchs, H.H. Woher, Nucl. Phys. A650 
(1999) 97. 

B. D. Serot, J. D. Walecka, Advances in Nuclear Physics, 
16, 1, eds. J. W. Negele, E. Vogt, (Plenum, N.Y., 1986) 
P. Ring, Prog. Part. Nucl. Phys. 78 (1996) 193. 

F. Weber, J. Phys. G25 (1999) R195. 

B. D. Serot and J. D. Walecka, Int. J. Mod. Phys. E 6 

(1997) 515. 

R.J. Furnstahl, B.D. Serot, and H.-B. Tang, Nucl. Phys. 
A615, 441 (1997); H.W. Hammer, R.J. Furnstahl, Nucl. 
Phys. A678, 272 (2000). 

M. Lutz, B. Friman, Ch. Appel, Phys. Lett. B474 (2000) 
7. 

G. J. Horowitz, B.D. Serot, Phys. Lett. B137 (1984) 287; 
Nucl. Phys. A464 (1987) 613. 

W. Botermans, R. Malfliet, Phys. Reports 198 (1990) 
115. 

B. ter Haar, R. Malfliet, Phys. Reports 149 (1987) 207. 

C. Fuchs, T. Waindzoch, A. Faessler, D.S. Kosov, Phys. 
Rev. C58 (1998) 2022. 

T. Gross-Boelting, C. Fuchs, A. Faessler, Nucl. Phys. 
A648 (1999) 105. 

C. Fuchs, H. Lenske, H.H. Woher, Phys. Rev. C52 (1995) 
3043. 

R. Machleidt, The Meson Theories of Nuclear Forces 
and Nuclear Structure, Advances in Nuclear Physics 19 
(1989) 189. 



R. Machleidt, |nucbth/991105c . 

C. Fuchs, E. Lehmann, L. Sehn, F. Scholz, J. Zipprich, 
T. Kubo, and Amand Faessler, Nucl. Phys. A603 (1996) 
471. 

N. Bastid (FOPI collaboration), proceedings to Structure 
of the Nucleus at the Dawn of the Century, (Publisher, 
Bologna 2000). 

C. Fuchs, P. Essler, T. Gaitanos, H.H. Woher, Nucl. 
Phys. A626 (1997) 987. 

R.K. Puri, N. Ohtsuka, E. Lehmann, A. Faessler, H.M. 
Martin, D.T. Khoa, G. Batko, and S.W. Huang, Nucl. 
Phys. A536 (1992) 201. 

L. Sehn, H.H. Wolter, Nucl. Phys. A601 (1996) 473. 

D. T. Khoa, N. Ohtsuka, M.A. Matin, A. Faessler, 
S.W. Huang, E. Lehmann, R.K. Puri, Nucl. Phys. A548 
(1992) 102. 

C. Fuchs, Amand Faessler, E. Zabrodin, Y.M. Zheng, 
Phys. Rev. Lett. 86 (2001) 1974. 



14 



[37] N. Ohtsuka, R. Linden, A. Faessler, F.B. Malik, 
Nucl.Phys. A465 (1987) 550; I. Izumoto, S. Krewald, 
A. Faessler, Nucl.Phys. A341 (1980) 319; 

[38] M. Baldo, I Bombaci, G. Giansiracusa, and U. Lombardo, 
Phys. Rev. C40 (1989) R491. 

[39] W. Zuo, I. Bombaci, U. Lombardo, Phys. Rev. C60 

(1999) 024605. 

[40] F. de Jong and H. Lenske, Phys. Rev. C54 (1996) 1488. 
[41] G.Q. Li and R. Machleidt, Phys. Rev. C48 (1993) 2707. 
[42] H. Sorge, H. Stocker, W. Greiner, Ann. Phys. 192 (1989) 
266. 

[43] S. Huber and J.Aichelin, Nucl. Phys. A573 (1994) 587. 
[44] Ch. Hartnack, R.K. Purl, J. Aichelin, J. Konopka, S.A. 

Bass, H. Stocker, W. Greiner, Eur. Phys. J. Al (1998) 

151. 

[45] V.S. Uma Mahoswari, C. Fuchs, Amand Faessler, L. 

Sehn, D. Kosov, Z. Wang, Nucl. Phys. A628 (1998) 669. 
[46] C. Fuchs, H.H. Wolter, Nucl. Phys. A589 (1995) 732. 
[47] A. Andronic et al. (FOPI collaboration), Nucl. Phys. 

A661 (1999) 333c. 
[48] C. Pinkenburg et al. , Phys. Rev. Lett. 83 (1999) 1295; 

P. Danielewicz, Roy A. Lacey, et al. , Phys. Rev. Lett. 

81 (1998) 2438. 
[49] V. Ramillien et al. (FOPI-Collaboration) , Nucl. Phys. 

A587 (1995) 802. 
[50] T. Gaitanos, H.H. Wolter, C. Fuchs, Phys. Lett. B478 

(2000) 79. 

[51] P. Crochet (FOPI-Collaboration), XXXIV International 
Winter Meeting of Nuclear Physics (I. lori, Bormio 
1996); N. Bastid et al. (FOPI-Collaboration), Nucl. 
Phys. A622 (1997) 573. 

[52] W. Reisdorf, Nucl. Phys. A612 (1997) 493. 

[53] F. Rami et al. (FOPI Collaboration), Nucl. Phys. A646 
(1999) 367. 

[54] E. Zabrodin, C. Fuchs, L.V. Bravina,. A. Faessler, Phys. 

Rev. C63 (2001) 034902. 
[55] A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C58 

1671 (1998) 1671. 
[56] A. Andronic et al. (FOPI-Collaboration), Phys. Rev. 

C64 (2001) 041604. 
[57] P. Crochet et al. (FOPI-Collaboration), Nucl. Phys. 

A624 (1997) 755. 
[58] T. Odch et al. (ALADIN Collaboration), Phys. Rev. Lett. 

84 (2000) 4557. 
[59] C. Sturm et al. (KaoS Collaboration), Phys. Rev. Lett. 

86 (2001) 39. 

[60] S. Typel, T. v. Chossy, W. Stocker, H.H. Wolter, pro- 
ceedings to Int. Symp. Nucl. Structure Physics, od. R. 
Castcn et al. (World Scientific, 2001). 

[61] F. Baffin, K. HagUn, and W. Bauer, Phys. Rev. C54 
(1996) 1375. 



15 



